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Abstract. We investigate distributions of hyperbolic Bessel processes. We 
find links between the hyperbolic cosinus of the hyperbolic Bessel processes and 
the functionals of geometric Brownian motion. We present an explicit formula 
of Laplace transform of hyperbolic cosinus of hyperbolic Bessel processes and 
some interesting different probabilistic representations of this Laplace trans- 
form. We express the one-dimensional distribution of hyperbolic Bessel process 
in terms of other, known and independent processes. We present some applica- 
tions including a new proof of Bougerol's identity and it's generalization. We 
characterize the distribution of the process being hyperbolic sinus of hyperbolic 
Bessel processes. 

1. Introduction 

The important role which functionals of Brownian motion play in many fields of 
mathematics (as mathematical finance, diffusion processes in random environment, 
probabilistic studies related to hyperbolic spaces etc.) is motivation to study the 
wide class of different diffusion processes connected to those functionals somehow 
(see for example [7], [5], [9], [14] and Q~5] ). In this work we consider the special 
class of such diffusions: Brownian motion with a special stochastic drift. It is 
well known in literature that processes like Ornstein-Uhlenbeck, Bessel processes or 
radial Ornstein-Uhlenbeck are examples of Brownian motion with stochastic drift 
connected to special functions: parabolic cylinder, Bessel or Kummer functions (see 
[3]). In this work we consider another interesting class of such diffusions - hyper- 
bolic Bessel process which are connected to Legendre functions. We establish the 
distribution of a hyperbolic Bessel process. Borodin presented the computation of 
the transition density function for hyperbolic Bessel process (see formulas (5.4) and 
(5.5) in [3]). The method used by Borodin relied on the connection between the 
Laplace transform of transition density function and two increasing and decreasing 
solutions of some ordinary differential equation of the second order (for the more de- 
tails of this method see also [5]). Gruet established the transition density functions 
of hyperbolic processes by the planar geometry method (see [11], [10]). However, 
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the results obtained by both authors are technically very complicated. The role of 
the hyperbolic Bessel processes in the world of Brownian motion functionals was 
mentioned also by Byczkowski, Malecki and Ryznar (see [5]). 

Our research of hyperbolic Bessel process is completely different (when compared 
to the mentioned) and purely probabilistic. We present results for hyperbolic Bessel 
processes with the index a > —1/2 and for the fixed time. We find the link be- 
tween a hyperbolic cosinus of hyperbolic Bessel process and functionals of geomet- 
ric Brownian motion. We present different probabilistic representations of Laplace 
transform of a hyperbolic cosinus of hyperbolic Bessel process R. We also give prob- 
abilistic representations of density of coshi? t . The link between hyperbolic Bessel 
process and functionals of geometric Brownian motion enable us to establish simple 
explicit form of Laplace transform of the vector (e Bt+kt , J e 2 ^ B " +ku ' , du) for a stan- 
dard Brownian motion B and an nonnegative integer k. These results are obtained 
for fixed t (not for time being stochastic) and that they can be effectively used in 
numerical computations. We express the distribution of hyperbolic Bessel process 
R in terms of squared Bessel process X and vector {B^\ A^) of Brownian motion 
with drift and the integral of geometric Brownian motion, which is independent of 
X. The joint distribution of (B^\ A^) is known in literature and is expressed 
via the density function (e.g. see [H]). As an interesting example of applications 
we find a new simple proof of the Bougerol's identity. We also characterize the 
distribution of sinh(i?) for R being hyperbolic Bessel process. It is interesting that 
the distribution of such process is deduced from the stochastic process which is in 
a sense a generalization of a squared Bessel process. At the end, we outline the 
case of time being exponential random variable independent of a Brownian motion 
driving the hyperbolic Bessel process. In this case many explicit and computable 
results are known in literature (see for example [5], |14| or |15|). 

2. Hyperbolic Bessel processes with fixed time 

We consider a complete probability space (f2, J- ', F) with filtration F = (-^)te[o,oo) 
satisfying the usual conditions and T = Too . We define a hyperbolic Bessel process 
R with index a G [—1; oo) as a diffusion, starting from a nonnegative x, given by 

1 r* 

(1) R t = x + B t + (a + -) / coth R u du. 

2 Jo 

There exits a unique solution to ((T|) defined up to explosion time and behavior of 
R depends on the index a. For a =/= —1/2, the point is an isolated singular point 
of (fl}. For a > — i, a solution of (TT]) is strictly positive for x > and is unique 
up to T a = inf{£ > : X t = a} for every a > x > (see Theorem 2.16]). 
In particular, for a > — ^ and x > there exits a unique solution defined up to 
the explosion time. For a = — i, a process R satisfying (fT]) is a Brownian motion 
starting from x. For a G [—1; —5), a solution of ([I]) is strictly positive for x > 
and negative for x — 0. It is unique up to T a = inf{< > : X t — a} for every 
a > x > (see [5J Theorem 2.17]). In particular, for a G [—1; — |) and x > there 
exits a unique positive solution defined up to the explosion time. 

Recall that Gruet ([TT]) defines a hyperbolic Bessel process with index a > —1/2 
as a nonnegative diffusion R with generator 

1 d 2 , 1 N , . . d 

(2) ^-2^ + ( a+ 2) COth( ^ 
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In the case of— 1/2 < a < the definition is completed by requirement that 
is an instantaneously reflecting point. It is also assumed that a starting point is 
nonnegative. The similar definition of hyperbolic Bessel process can be found in 
Borodin [2]. In [TT], by the identification of the Green function, it is stated that 
a reasonable candidate for hyperbolic Bessel process with a = —1/2 is a reflected 
Brownian motion (see the end of Section 3) . We define a hyperbolic Bessel process 
with a = — 1/2 as a Brownian motion instead of the reflected one, since a Brownian 
motion is a solution of the (flj . Such an extension of definition of hyperbolic Bessel 
process by a Brownian motion allows to recover the interesting results about the 
distribution of J*e 2Bu du (including the Bougerol's identity). 

To start the discussion about the properties of hyperbolic Bessel process we 
define a process 9 being the solution of the SDE 

(3) d0 t = yJ\0?-l\dB t + (a + l)\O t \dt, 

where a > — 1, and such that 6q = x > 1. Observe that the diffusion coefficient 
a(x) = \/\x 2 — 1| is locally Lipschitz. Moreover, for b(x) — (a + \)\x\ and for x € M, 

\a(x)\ + \b(x)\<(2 + a)(\x\+l). 

Thus, the SDE (J3J has an unique strong non-exploding solution (see [5J Chapter 
111,2])). Now, we consider the diffusion r\ t := 9 t — 1. It is a diffusion with drift 
and diffusion coefficients equal to b(y) = (a + l)\y + 1| and a(y) = y/\y 2 + 2y\, 
respectively. We observe that the point is an isolated singular point for r\. For 
a > and a > we have 

Texpf r 2 Mdy]dx= r(^Y +1 dx>a^ f x~^dx = 00. 
Jo a 2 (y) > Jo \x 2 + 2x) ~ Jo 

Therefore, from Theorems [51 Theorem 2.16] and [5J Theorem 2.17] follows that 
r]t > for all t > 0, provided 770 > 0. So, if 9o > 1, then 9 t > 1 for all t > 0. Hence, 
the process 9 satisfies for a > the SDE 

(4) d9 t = ^9 2 - 1 dB t + {a + l)9 t dt. 

It turns out, that in the case of a > we can recover a hyperbolic Bessel process 
from the process 9. As usual, by ar cosh we denote the inverse function of cosh. 

Theorem 2.1. If a > and x > 1 the process Rt — arcosh(9t) is a hyperbolic 
Bessel process, where 9 is the solution of 

Proof. Observe that, by definition, 

(5) Rt = ar cosh(0 t ) = In (o t + ^9 2 - l) > 0. 

As 6> t > 1 for alH > 0, we can use the Ito lemma to obtain that R satisfies (JXJ) . □ 

Proposition 2.2. A hyperbolic Bessel process R with index a > —1 does not 
explode. 

Proof. For a > we use Theorem 12.11 which implies that R does not explode 
since 9 does not explode. For a = — | the process i? is Brownian motion. For 
a E [—1; 0) \ { — 5}, we observe that cosh(i?) and 9 have the same generator Ag on 
C 2 , the space of twice continuously differentiable functions with compact support. 
From the uniqueness of solution of martingale problem induced by Ag we conclude 
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that cosh(i?) and 9 have the same distributions. The diffusion 9 does not explode, 
and consequently R does not explode. □ 

Now we define a process T^ a \ which appears in the sequel. Let B be a standard 
Brownian motion, a G R, A > and 

B t +at 

(6) 



r(«) - 
1 1 — 



1 + A / e B ^+ au du 



There are connections between processes r'~5) and cosh of a hyperbolic Bessel 
process with a = — 1. Namely, in Jakubowski and Wisniewolski P21 Theorem 2.4] 
it is proved that for a hyperbolic Bessel process R with index a = —1, Ro = x > 0, 
and A > we have 



(7) 



£ e -Acosh(fl t ) _ e -Ajg e -A(cosh(a;)-l)r t 



Now we present a form of Laplace transform of cosh of hyperbolic Bessel processes 
with a > —1/2, one of the most important results of this paper. 



Theorem 2.3. Let a > —1/2 and t > be fixed. If R is a hyperbolic Bessel process 
of the form §Q with x > 0, then for A > 

/ A 2 rt 
exp y — Acosh(a;)V t — 



(8) 



3 



exp ( — A cosh R t 



V„fdu 



where V t — e \ a+ 2 / t+Bt and B is a standard Brownian motion. 
Proof. Set 9t = coshi? t > 0. So, 9q = cosh a; and 

(9) d9 t = \j9 2 t -l dB t + aJOtdt, 

where a = a + 1. From the Ito lemma we obtain 

(10) 



I \2 

de- X9t = -\e- xet ^9 2 t ~ UB t - Xae~ X9t 9 t dt + —e~ X8t (9 2 - l)dt. 



Now, we show that J Q e xe " \j9\-\ dB u is a martingale. It suffices to prove that 



(11) 



E 



;du < oo 



for every t > 0. From © we have 

c-l 



'o 



< 3 



y/61 - ldB u + a I 9 u du 
t 

y/ei - ldB u 









J 9 u duj 







Using estimation of stochastic integrals and the Jensen inequality for f(x) = x 2 we 
deduce that 



E 









) 2 +« 2 ( 


j o u d U y 


< CE f 






Jo 



t u - l)du + a 2 ti 



9ldu. 



Taking C\ = 30 2 , C% = 3(C + ai) with a\ > a 2 T, we conclude that 

E9 2 <d+ C 2 E / 9ldu, 
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for t £ [0, T], so by the Gronwall lemma 

E6> t 2 < Cie C2t , 

and (fTTj) follows. 

Therefore, from (flU)) , we infer that 



2 

A6>„ 



Ee- xet = e- xe ° - aX f E[e' xe -e u }du+^- f E[e' xo " (0 2 u - l)]dit 



Define A) := Ee X6t , so the function p is bounded. Using p"2")) . we deduce that 
p satisfies the following partial differential equation 

dp A 2 dp A 2 <9 2 p 

(13) d-t=-T p + aX d\ + Td\*> 

withp(0,A) = e- xe °. Set V t = V e( a+ ^ t+Bt . Then 

dV t = V t dB t + aVtdt, 
with a = a + 1 . The generator of V is of the form 



2 J2 



V 2 da; 2 dx 



Observe, by {EJ), that p € C^ 2 ((0, oo) x (0,oo)) and from the Feynman-Kac theo- 
rem (see |13[ Chapter 5 Theorem 7.6]) we know that the only bounded solution of 
partial differential equation 

^ = Avu-^-u, u(0,x) = e- e ° x , 
i.e. (fT3|) . admits the stochastic representation 

u(t, x) = Eexp ( - 9 V t - i J V*du) , 
where Vb = x. This implies assertion of the theorem. □ 

As an easy consequence we obtain 
Proposition 2.4. The Laplace transform of the vector (e Bt , J Q e 2B "du) is given by 

(14) Eexp(-7e St - y J e 2S "du) = Eexp ( - Acosh(a; + B t fj , 



/or 7 > and A > 0, where x = ar cosh j. 

Proof. Proposition follows by applying Theorem 12.31 with a = — g- D 

Remark 2.5. The form of the Laplace'a transform of the vector (e Bt , J Q e 2Bu du) 
in the last proposition, i.e. (|14p . is very simple when it is compared to the form 
of density of this vector obtained by Matsumoto and Yor [T3] . Indeed, the density 
given in |14| has the oscillating nature in the neighbourhood of and is not conve- 
nient for computational use (see [3]). The knowledge of Laplace transform make it 
possible to invert it numerically and obtain the density of (e Bt , J e 2B "du). It is also 
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important to realize that it enables to obtain the density for (e s ' M) , J* e 25 " 1 ' du), 
where := B t + fit, \l € R, as we know from P3] that 

¥(b^ ] G dx, e 2B ^du G dyj = e" x -^ t/2 ¥^B t G dx, J e 2B "du G dy 

Theorem 2.6. Let k G N, 7 > 0, A > and R be a hyperbolic Bessel process of 
the form 

(15) Rt — ar cosh(7/A) + B t + k / cothi? u du. 

Jo 

For A > we have 



E 

where 



exp(-Acoshi ?t )=(-l)^^l 



p(j,X) = Ee~ Acosh ( a ' rcosh ( 7 / A )+ s *) i 
Proof. By Proposition 12. 4\ for a; = arcosh( 7 /A), 

p(7, A) = Ee- Acosh ^+ s *) = Eexp ( - 7 e B ' - y jf* e 2B "dn 
Define the new probability measure Q by 

dQ = 

dP T t 

Then Vt = Bt — kt is a standard Brownian motion under Q and 

/ fc 2 A 2 /"* \ 

p(7, A) = Eq exp ( - kB t + —t - 7 e Bt - — J e 2B »du) 

= e**E Q exp ( - k(V t + kt) - 7e W+ fc ') - y T e 2(y "+ fctl) dn) . 



The result follows from Theorem 12.31 after taking fc-th derivative of p with respect 
to 7. □ 

Proposition 2.7. For k G N, A > 0,7 > and Y~ t = e B t+ w 

* /"Wa_, ^„ 4*^(7, A) 



Eexp ( - 7 y t - y jf K u 2 d u ) = (-l) fc e^ 



c>7 fc 



w/iere 

p( 7 ,A) =Ee- Acosll(:E+i3t) 

and x = ar cosh(j). 

Proof. Proposition follows from Theorems 12.61 and 12.31 □ 



Theorem 2.8. If R is a hyperbolic Bessel process of the form ([I]) with a > —1/2, 
then for A > 

Eexp(- Acoshi? t ) = ^e-^E(l {Vt > lBtn V t h t {B t ,x)J a (\0(B tl V t ))y 

where a = a + ^, Jq is the Bessel function of the first kind of order 0, h t (z,x) = 
e £+az-\cosh(x)e^ ^ z ^ = y /2 e x / 2 (cosh z - coshx) 1 / 2 , z > \x\ and V, B are two 
independent standard Brownian motions. 
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Proof. Let Wt = Bt + at and let the new probability measure Q be given by 

2 

-oBi-T-i 



From Theorem I2J3I we have 



Eexp ( - Acoshi? t ) = Ecxp ( — A cosh(x)e a * +St - y / e 2[Bu+au) du 

( a 2 A 2 "* 
= Ejjexp [aWt — —t — Acosh(ai)e Wt — 

2 / A 2 /*' \ 

= e^'Eexp [aB t - Acosh(.x)e St - — J e 2Bu duj, 



2W -du 



Eexp (-Acosh R t ) = e"^*E e^-Acosh^'jgYg-^ /* e SB «<i«| B( j 
To finish the proof we use the conditional Laplace transform (see (5.5) in |14j ) 



E(e -2 2-^ e2B " du |-B 4 = sc 



1 6 2t 



e 2' Jo(X(f>(x, z))dz. 



□ 



Theorem 2.9. If R is a hyperbolic Bessel process of the form ([T]) with a > —1/2, 
then for A > 

Eexp ( - Acoshfl,) = e -A E [ e -A(co 8 h(x)-i)if +4) (l + \J* e^+^+^du)^ 

where B is a standard Brownian motion and f[ is given by (J5J) • 
Proof. Set Y^ = e Bt+at . Define the new measure Q by 



Theorem 12.31 implies 



Ee -Acoshi? f = jg 



Er 



= Er 



- cosh(i) AY t <Q+ 2 '+A /„' y!° + 2 ! dB„ 



cosh(a:)AY t < " +7) +A(Y t ( ° +5 '-l)-(a+l)A/ t 



DA f„ t K (Q + '' 



where we have used the fact that 



(16) 



Y, 



(a+i) 



1 



(«+§) 



Yu dB u + (a + 1) / r u v " T2; du. 



n 



From the Girsanov theorem, the process 14 = _B t + A Ki Q+2 ''c?u is a standard 
Brownian motion under Q. By the result of Alili, Matsumoto, and Shiraishi PJ 
Lemma 3.1], we have 

(17) 



B t = V t -hi(l + xJ e v " +{a+ ^ u du 
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Thus 



and 



(«+*) 



,y t + (a+i)t 



I = e 



-At, 



-A(cosh(a:) — 1)- 



-(a+l)(V t -B t ) 



,-At, 



-a-l 



□ 



Considering a = —1/2 in the previous result we obtain 
Corollary 2.10. For A > 0, x > 0, 



E 



exp ( — A cosh(_B t + a;) 



e E 



e-ACcosh^-i^Yi + A / e B -du 



-1/2 



In the sequel we use the following notation: 



(18) 



.4 



(n) 



5 2(B U+Q u) , a 



with £> being a standard Brownian motion. 



Theorem 2.11. If R is a hyperbolic Bessel process of the form ([!} with a = and 
x = 0, i/ien i/ie density o/coshi? f on [l,oo), t > 0, has the form 



(19) 



P(coshi?i £dz)=E 



1 



4A 



(i) 



exp 



z - 1 



<iz. 



Proof. Fix i > 0. Using Theorem 12.91 for a = x = and A > we obtain 

1 

1 + A J*e B ^du- 

From the last equality and the scaling property of Brownian motion we have 



g e -A(cosh_R t -l) _ jg 



(20) Ee -A(coshB 4t -l) =JE 



1 



-1 + 4A J c ]e B ^+ 2u du 
From Theorem 2.8 in [T2] we know that 



= E 



1 



1 + 4AA 



(i) 



E 



1 




= 1 - 4A f 


1-1 + 4A^ 1} - 


Jo 



where Gt(y) = Ele 



Since G t is differentiable, integration by parts formula yields 
1 



E 



1 + 4AA 



(i) 



1+/ G t {y)(e-^y)'dy 



G>(y)e- 4X ydy. 
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From the last equality and the definition of Gt we have 



E 



1 



1 + AXA 



(i) 



-4Ay E 



1 



.4 



exp 



(21) 



e _Az E 



4A 



(i) 



exp 



4A 



(i) 



From {2D) and (J3TJ) we conclude that 
(coshi?4 t — 1) G dz 

for z > 0, which ends the proof. 



4A 



(i) 



exp 



4A 



(i) 



<iz. 

dz, 



□ 



Remark 2.12. Using the explicit form of the density of Ap and Theorem I2.11[ 
we can obtain the integral form of density of coshi? t , with a = and x = 0. 

Proposition 2.13. If R is a hyperbolic Bessel process of the form (TTJ) with a = 
and a; = then on [1, oo) 



P(coshi? t G dz) 



H(¥) 



dz, 



G t (y) = e-*/ 2 Eexp (i? t + ^(s t 2 - ^ 2 (5 t ) 



(22) 

/or t > 0, where 
(23) 

B is a standard Brownian motion and 

(24) </ 3 y( z ) = hi (?/ e ~ Z + cosh(z) + \J y 2 e~ 2z + sinh 2 (z) + 2ye~ z cosh(z)^ . 

. y v 

Proo/. As before, let G t (t/) = E e ^ . As an easy consequence of the Matsumoto- 

Yor result [Ffl Theorem 5.6] we obtain that Gt is given by (j2"5)) (see the proof of 
Theorem 2.5 in Q~2]). Proposition follows from Theorem 12.111 and observation that 



E 



AA 



(i) 



exp 



z - 1 



□ 



(law) 



In the next theorem we prove that Rt - l — ' F(B, A, X) for some functional F, so 
the distribution of R t can be presented as functional of A, B and a squared Bessel 
process X of index a independent of B. 

Theorem 2.14. If R is a hyperbolic Bessel process of the form {l} with a > —1/2, 
Rq = x, then 

(25) cosh( J R t ) (1 = ) l + 2A^ +1) A 1 

for every t > 0, where X is a squared Bessel process of index a, independent of a 
standard Brownian motion B and starting from i(cosh(a;) — l)e 2B */ 4 /A^" +1 \ 
Moreover, for a > —1/2, 

(law) 



(26) 



Rt 



1 + 2Af / a A +l) X l 



10 



JACEK JAKUBOWSKI AND MACIEJ WISNIEWOLSKI 



Proof. From Theorem 12.91 and the form of Laplace transform of squared Bessel 
process (see \17\ Chapter XI, p. 441]) we have, for arbitrary A > 0, 

Ecxp ( - Acoshi? t ) = e-A^-AtcoshM-Drr^ + c B ^ a+ ^ u du) 

JliCXp^ AC A f(A+ (a + l)»,, i/( ,S t + (a + l)t, ! 



e- x 



where X is a squared Bessel process of index a, independent of a standard Brownian 
motion B and starting from (cosh(x) — 1). It is now clear that 

(27) cosh(i? t ) ( '= > 1 + e*^^^^.^^^. 

Using the scaling property of squared Bessel process (see [17j Chapter XI, Prop. 
1.6]) we obtain 

(28) 1 + e^^X ro eBu+(a+ ^ d „ /(2eBt+(Q+i)t) =l + ^i { e^^du, 

where AT is a squared Bessel process of index a independent of B and starting 
from the stochastic point X = 2(cosh(a:)-l)e Bt+ ( Q+ ^7 f* e B ^ +( - a+ ^ u du. Hence, 
taking At instead of t, we infer that 

cosh(i? 4t ) (1 = } 1 + \lCx e B ^ a+ ^du (1 = } 1 + 23f a+1) X u 
2 Jo 

where X is a squared Bessel process of index a independent of a standard Brownian 

— — , . , . ... , ,n(2o+i) ,— r (2a+l) , —r(2a+l) 

motion B and starting from Xq = i(cosh(.x) — l)e 4 , where A 4 

is given by (fT8]) with i? instead od B. The proof is complete. □ 

Remark 2.15. Fix a > —1/2. For every t > 0, from the proof of Theorem 12.141 
we obtain that 

(29) Rt ( '= } ar cosh (l + {l/2)a t X^\ , 

where a t — J Q e Bu+i - a+ ^ u du, X is a squared Bessel process of index a > —1/2 
independent of a standard Brownian motion B and starting from the random point 
X = 2(cosh(x) - l) e St +( Q +5)7 J* e B "+^ a+ ^ u du (see $7$ and (ggj)). 

The last theorem, in the special case x = 0, gives 

Proposition 2.16. If x = and a > —1/2, £/ien /or every t the density function 
of cosh(i? t ) on [l,oo) /ias the form 

(30) P(cosh(^)e^) = i i T ^^E 



-( Z -l)/(4A< 2 / ° + 1) ) 



a+l 



(7.-J . 



If x = and a > —1/2, £/ien for every t the density function of R t on [0,oo) has 
the form 
(31) 

m e <fc) = JL (cosli(7 -l) a sinh(.) 
v ; 4"+! r(a + l) 



-(co S h( z )-l)/(4A< 2 / ° + 1 ») \_ 



a+l 
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Proof. From Theorem 12. 141 for x = and a > — 1/2 we obtain 

(law) 



(32) 



cosh(i? t ) 



where X is a squared Bessel process of index a, independent of B and starting from 
0. For a > -1/2 it holds that 

(law) 



(33) R t 
Hence for x = and a > —1/2 

P(cosh(i? i ) < z 

and for .t = and a > —1/2 



1 + 2Af / a A +l) X 1 



X 1 < 



z-1 



2A 



(2q+1) 
t/4 



\Rt <z)=F[X 1 < 



cosh(z) — 1 

(2a+l) 
t/4 



2A 



From the last equality and the form of the density of squared Bessel process (see 
(TTI Chapter XI, Cor. 1.4]) we have for x = and a > -1/2 



I (cosh(i? t ) G dz) 



l\«+i (z - l) c 



T(a + 1) 



E 



3 _ (z _l)/(4Af / :+ 1 ')_ 



1 



.4 



;2o+i; 

t/4 



a+1 



dz. 



and for x = and a > —1/2 
P(i? t G dz) = 

■ 1 ^ "+ 1 (cosh(z) - l) a sinh(z) 
r(a + 1) 



S 



s -(cosh( Z )-l)/A^ +1 >_ 



a+l 



dz. 



□ 



Remark 2.17. Using the explicit form of the density of and Proposition 

12.161 we can obtain the integral form of density of R t for a > —1/2 and Rq = 0. 
This formula is new and differs significantly from the density obtained by Borodin 
(see formula (5.4) in [4]). 

The next two facts follows immediately from 12. 141 and in slightly different forms 
can be found in [T4"] . 

Corollary 2.18. For any z > 1 we have 



E 



1 



exp 



(cosh(z) — 1) 



2 v /cosh(z) - 1 



t sinh(z) 

Proof. It follows from Proposition 12.161 for x — and a — — ^ and the form of 
density of cosh (St). □ 

Proposition 2.19. For every A > we have 

A _ E(cosh(i?4t) - 1) A 
* 2 A E(£i) A 
where B is a standard Brownian motion. 

Proof. It follows from Theorem 12.141 for x = and a = — 1/2. □ 
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(35) sinh(.T + B t ) {l = ] sinh(a;)e iJt + W Af ■ 



In the next proposition we present an alternative purely probabilistic proof of 
generalization of Bougerol's identity (see [T|, [21 Prop. 6.1] or [M) Section 3]). 

Proposition 2.20. Let B and W be two independent standard Brownian motions. 

a) [Bougerol's identity] For t > 0, 

(34) sinh(i? f ) ( = W) W^ t . 

bjlfxeR, then 

(35) 
Proof, a) 

(W A y °= > A t { Wl f °= W) C ° sHB f - 1 = e2kBit + \ 2kBit 2 (1 = } (sinh(5 t )) 2 , 

where in the second equality we use Theorem 12.141 with x = and a = —1/2. 
Hence, using the fact that Wa ± and sinh(B t ) are symmetric random variable we 
obtain ([M]) . 

b) By the same arguments as in a), we conclude that 

| sinh(a; + B t )| (1 = } | sinh(a;)e Bt + W At \ 
and we can expect that we can omit modulus. Indeed, define 

M t = e St (sinh(x) + / e - B "dW u ). 
Jo 

Observe, by Proposition 2.1 in [9], that 

sinh(a;)e Bt + W At (1 = ) sinh(a;)e St + e Bt [ e~ B "dW u = M t . 

Jo 

It is easy to check that M is a diffusion with the generator 

1 , 9 n d 2 1 d 

It is also evident that SDE corresponding to the above generator has a unique 
strong solution. To finish the proof it is enough to observe that Mt := sinh(x + Bt) 

is a diffusion such that generators of M and M are equal on C%, so M M. □ 

In the next proposition, we deduce from Theorem 12.201 the simple form of char- 
acteristic function of vector (e Bt , W At )- 

Proposition 2.21. Let uel. Then, for v ^ 0, 

-^ e iue Bt +ivWA t _ ggTOsinh(arsinh(u/ii)+J3 t ) 

Proof. From Theorem 12.201 we infer that 

-e Bt + W At ( '= ) sinh(ar sinh(u/u) + B t ), 
v 

which implies (|3T)) . □ 

Now we establish the representation of law of the process sinh(i?), where R is a 
hyperbolic Bessel process with index a = 0, in terms of functionals of independent 
Brownian motions. 
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Theorem 2.22. Let x > and R be hyperbolic Bessel process starting from x and 
with the index a — 0. Then 

(smh(Rt),t> 0) (I = } 

(38) (e- B ' +t / 2 (sinh(a;) + j\ B ^ 2 dV u ) 2 + ( J* e B -- u ' 2 dZ u ) *) ' , t > o), 

where B, V, Z are three independent standard Brownian motions. 

Proof. Define £t = sinh 2 (i? t ). Then the diffusion £ satisfies the following SDE 



and the generator of £ is of the form 
(39) At = 2(x 2 +x)-^2 + (2 + 3x)- rl 



dx 2 dx 

Let's denote the diffusion and drift coefficients of £ by a(x) = 2v 'x 2 + x and (i(x) — 
2 + 3a:. Observe that for any x > and y € (x — 1, x + 1) 

(a(x) - tr(y)) - 4(^ 2 + ic - + y) 2 < 4(|x 2 - y 2 \ + \x - y|) 

= 4|x - y|(|a5 + y| + 1) < 4|x - y|(2|a?| + 2) 

< 4|a: - y|(o-(a;) + 2) < 4jac - ?/|(cr 2 (x)/2 + 5/2). 

The uniqueness of the solution for the corresponding SDE follows now from |17[ 
Chapter IX Ex. 3.14]. Thus, there is a unique solution of a martingale problem 
induced by A%. 

For a standard Brownian motion B we define, as previously, Y t = e Bt ~ t l 2 . We 
consider now a SDE of the form 

(40) dX t = 2yfx t Y t dW t + 2Y 2 dt, 

where Xq = sinh (x) and W is a standard Brownian motion independent of B. 
Observe, using the ltd lemma, that the process ip = yr is a diffusion with the 
generator A^p having the same form as A^ on on C 2 . Hence, by uniqueness of 
solution of the martingale problem induced by .4^, we obtain that the processes 
|sinh(i?)| and \y/X/Y\ have the same law. To skip the absolute value we realize 
that both processes are nonnegative. Moreover, the SDE (|4"0)) has a unique weak 
solution. Indeed, using the change of time r t = inf{.s > : J Q Y 2 du > t} we obtain 
that the process X T is a square of 2-dimensional Bessel process, and in consequence 
it is not difficult to see that the process X inherits "good" properties of a square of 
2-dimensional Bessel process. To finish the proof we observe that the unique weak 
solution of (|40p can be written as 



(41) X t = (sinh(i) + / YudVuY + ( / Y u dZ u y. 

Jo Jo 

To verify that (|4ip satisfies (|40[) it is enough to realize that W defined as 
r* (sinh(x) + J Q S Y u dV u )dV s + (/ S Y u dZ u )dZ s 

1/A , 



(42) W t := / l {As7t0} - 

Jo ^/(smh(x) + /; Y u dV u ) 2 + (J S Y u dZ u ) 2 

where A s = (sinh(a;) + f° Y u dV u ) 2 + (J° Y u dZ u ) 2 , is a standard Brownian motion. 
Simple use of the ltd lemma finishes the proof. □ 
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It turns out that the methods of the proof of Theorem 12.221 allows also to find a 
representation of law of hyperbolic Bessel process with index a > —1/2. 

Theorem 2.23. Let x > and R be hyperbolic Bessel process starting from x with 
the index a > —1/2. Then 



(43) 



'smh{Rt),t > 0) ( = w) (yX t /Y t ,t > o), 
where Y t = e s t -(«+i/2)« ; Xf satis fi es t h e $DE 
(44) dX t = 2^X t Y t dW t + 2(1 + a)Y?dt, 

Xq = sinh 2 (x) and £>, W are two independent standard Brownian motions. 
Proof. The proof goes in the same way as the proof of previous theorem. 
Theorem 2.24. Let X and Y be as in Theorem ^. 231 For any w > 0, t > 



(45) 



□ 



(^y 1 < w) = e[f s ( / Yldu,wY t 



where Fs(t,x) = ¥(St < w) and S is a Bessel process of dimension 2(1 + a) with 
So = sinh(x). 

Proof. Observe that the process X Tt is a square of 2(1 + a)-dimensional Bessel 
process for r t = inf{s > : J Q Y^du > t}. Observe also that conditionally, under 
the knowledge of trajectory of (Y s ,s < t), the process X Tt is still a square of 
2(1 + a)-dimensional Bessel process. Hence 



/Xt 
Yt 



E 



S f t Y 2 du < wY t 



o-{Y u ,u<t) 



( 



Yjdu,wY 



□ 



Proposition 2.25. Let k G N, x > and Y be a given continuous process such 
that E Jq* Y^du < oo for any t > 0. Then the unique solution of the SDE 



(46) 



dX t = 2JX t Y t dW t + kY t 2 dt, 



where W is a standard Brownian motion independent of Y and Xq = x, is of the 
form 



(47) 



t , 2 

YudWl 



E 

i=2 



Y u dWl] 



where W l , i — 1, . . . , k, are independent standard Brownian motions. 



Proof. As previously, uniqueness of solution of (|4^]) follows from the fact that the 
standard change of time (r t = inf{u : J Q Y 2 ds > t}) results in X becoming a square 
of k-dimensional Bessel process. We complete the proof by the direct checking 
that X defined by (|4T|) satisfies (|4"6")l using analogous arguments as in the proof of 
Theorem COS □ 
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2.1. Hyperbolic Bessel processes with stochastic time. 

Proposition 2.26. Let T$ be a random variable with exponential distribution with 
the parameter S > 0. Assume thatT\ is independent of a standard Brownian motion 
B. Then for a hyperbolic Bessel process of the form {1} with a > —1/2, we have 

pco roc 

(48) Ecxpf - Acoshi?T 5 ) = / / e- Xucosh{x) -? v p~<(u,l,y)dydu, 



where for 7 = ^25 + (a + , 

(49) pHu,1, V ) = -^e-^/^lJy) 

y-2 a u vu/ 

and I-y is modified Bessel function. 

Proof. We use Theorem 12.31 and the result of Matsumoto-Yor |14[ Theorem 4.11]. 

□ 
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